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ABSTRACT 

We present a careful analysis of the point source detection limit of the AKARI All- 
Sky Survey in the WIDE-S 90 fim band near the North Ecliptic Pole (NEP). Timeline 
Analysis is used to detect IRAS sources and then a conversion factor is derived to 
transform the peak timeline signal to the interpolated 90 ^m flux of a source. Combined 
with a robust noise measurement, the point source flux detection limit at S/N > 5 
for a single detector row is 1.1 ± 0.1 Jy which corresponds to a point source detection 
limit of the survey of ~0.4 Jy. 

Wavelet transform offers a multiscale representation of the Time Series Data 
(TSD). We calculate the continuous wavelet transform of the TSD and then search for 
significant wavelet coefficients considered as potential source detections. To discrimi- 
nate real sources from spurious or moving objects, only sources with confirmation are 
selected. In our multiscale analysis, IRAS sources selected above Aa can be identified 
as the only real sources at the Point Source Scales. We also investigate the correlation 
between the non-IRAS sources detected in Timeline Analysis and cirrus emission using 
wavelet transform and contour plots of wavelet power spectrum. It is shown that the 
non-IRAS sources are most likely to be caused by excessive noise over a large range 
of spatial scales rather than real extended structures such as cirrus clouds. 

Key words: methods: data analysis - surveys - infrared: galaxies - infrared: ISM. 



AKARI (previously known as ASTRO-F) is a Japanese in- 
frared astronomical mission launched in February 2006 (Mu- 
rakami et al. 2007). Its primary goal is to survey the entire 
sky with wide spectral coverage and high spatial resolution. 
It has two focal plane instruments, the Far-Infrared Surveyor 
(FIS, Kawada et al. 2007) covering a wavelength region from 
50 - 180 fim, and the Infrared Camera (IRC, Onaka et al. 



1 INTRODUCTION 



2007) covering a wavelength region from 1.8 - 26.5 fim. The 
FIS observes the sky in four photometric bands, referred to 
as N60, WIDE-S, WIDE-L and N160. In Tabled we give 
band name, band centre, detector array format and detector 
pixel size for each band. In comparison, IRAS has a pixel 
size of 1.5x4.7 arcmin^ at 60 ^m and 3.0x5.0 arcmin^ at 
100 fim. 



* E-mail: lingyu.wang@imperial.ac.uk 



The first phase of the AKARI All-Sky Survey started 
in May 2006, lasting for one half year. About 70% of the sky 
was observed by at least two independent scans during this 
period. In the subsequent complementary phase which ended 
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Table 1. Specifications of the Far-Infrared Surveyor (FIS). 

BAND N60 WIDE-S 



WIDE-L 



N160 



band centre 
detector array format 
detector pixel size 



65 

20x2 
26.8x26.8 



90 

20x3 
26.8x26.8 



140 

15x3 

44.2x44.2 



160 

15x2 

44.2x44.2 



[row X column] 



on August 26th, 2007, due to exhaustion of helium, the sky 
coverage was increased to about 94%. Similarly, IRAS sur- 
veyed 98% of the sky with a single scan, 96% with two con- 
firming scans and 72% with three or more scans (Moshir et 
al. 1992). 

The pre-flight 5(7 point source detection limit for a sin- 
gle scan (with 2 or 3 detector rows) was estimated to be 
0.6, 0.2, 0.4 and 0.8 Jy in the N60, WIDE-S, WIDE-L and 
N160 band respectively. The current in-orbit estimate based 
on noise measurement and absolute calibration using well- 
modelled objects is 2.4, 0.55, 1.4 and 6.3 Jy, with uncer- 
tainties 30% (Kawada et al. 2007). Thus, there is a fac- 
tor of 3 or more degradation in all bands caused by ex- 
cessive in-orbit noise, frequent glitches etc. On the other 
hand, the IRAS Point Source Catalogue (PSC) is complete 
to about 0.6 Jy at 60 /^m and about 1.0 Jy at 100 /xm in 
unconfused regions. At a 90% completeness, the IRAS Faint 
Source Catalogue (FSC) achieved a depth of ~0.2 Jy at 
60 jim. by coadding the data after point-source filtering. The 
IRAS 100 Atm band to which the AKARI WIDE-S 90 ^.m 
band is closer was not considered for catalogue qualification 
because of severe cirrus contamination. For the same reason, 
all IRAS flux densities at 100 pm were declared to be mod- 
erate quality. The FSC usually contains sources with 100 yLxn 
flux densities above 1 - 2 Jy. 

In this paper, we focus on confirmed detections from 
the timeline^ to give a reliable and updated point source 
detection limit estimate in the WIDE-S 90 ^m band, taking 
advantage of the improved data reduction pipeline. Section[2] 
describes the AKARI Performance Verification (PV) phase 
data near the North Ecliptic Pole (NEP) , robust noise mea- 
surement and 90 /im flux interpolation. Section [3] describes 
the Timeline Analysis processing, the conversion factor and 
the association between the seconds-confirmed (SCON) and 
hours-confirmed (HCONjEI sources which are not found in 
the IRAS catalogues and infrared cirrus emission. A wavelet 
multiscale analysis of the selected IRAS sources is presented 
in Section m which also includes an investigation on the ex- 
tended nature of the non-IRAS sources. Finally, we give our 
conclusions in Section (5) 




Figure 1. The area covered by the AKARI PV phase data near 
the NEP. 



from AprU 24th 2006 to May 7th 2006. The data file is in a 
dedicated FITS format called Time Series Data (TSD) for- 
mat. It has been processed and calibrated by the May 2007 
version of the Green Box which is part of a dedicated data 
reduction pipeline for the AKARI All-Sky Survey (Yama- 
mura et al. in prep.). The positional information is provided 
by the ground-based attitude determination system (GADS) 
and the rms positional accuracy is ~30". Eventually, the po- 
sitional accuracy will be improved to a few arcseconds using 
ESA's pointing reconstruction processing. The data covers 
approximately 40 square degrees (Fig. [l| and ~70 sources in 
the IRAS FSC with 60 /im fiux densities ranging from ~0.2 
to 10.1 Jy, excluding the Cat's Eye Nebula. This group of 
objects forms an important base to estimate the sensitivity, 
reliability and completeness of the AKARI All-Sky Survey. 

To compare with previous absolute calibration of the 
AKARI All-Sky Survey data, we have used another data 
set containing 9 asteroids whose positions, observed epochs, 
and predicted monochromatic 90 /im fiux densities based on 
model spectral energy distributions (SEDs; Mueller & Lager- 
ros 2002; Mueller et al. 2008) are listed in Table H 



2 THE DATA 

We have used 412 pieces of scan data near the NEP each of 
which is of one minute length, taken during the PV phase 

^ The one-dimensional data stream in a single detector pixel is 
defined as a timeline. 

^ Following the nomenclature of IRAS, seconds-confirmed sources 
are detected in at least two detector rows in a single scan; Hours- 
confirmed sources are detected in at least two independent scans 
which are hours apart from each other. Confirmation is essential 
to detect true fixed sources. 



2.1 Robust noise estimation 

We use a 3(t clipping method to estimate the Gaussian noise 
in the data. The main processing steps are: 

(1) For each detector data stream, we calculate and sub- 
tract the median with a specified width of 192 discrete sam- 
plings which corresponds to ~8 seconds or ~28.8 arcmin 
(background subtraction). 

(2) We discard fiagged/bad data which include signals 
caused by resetting, calibration lamp, dead pixels, cosmic- 
ray glitches, etc. In addition, we delete any region affected 
by the Cat's Eye Nebula which is a very bright extended 
source. 

(3) A 2-point boxcar smoothing is applied to the time- 
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Table 2. Asteroids used in absolute calibration. 



ID 


NAME 


RA (J2000) 


DEC {J2000) 




DATE 


TIME 


0210 


6 Hebe 


309.970 


-8.744 


13.197 


2006.04.30 


14:55:29 


0220 


511 Davida 


315.824 


-19.530 


11.373 


2006.05.03 


04:02:44 


0240 


1 Ceres 


328.007 


-21.064 


172.95 


2006.05.14 


00:51:00 


0295 


230 Athamantis 


202.658 


-11.580 


3.419 


2006.07.18 


01:56:30 


0480 


511 Davida 


307.679 


-28.328 


12.997 


2006.10.26 


16:20:22 


0510 


88 Thisbe 


133.914 


15.669 


7.271 


2006.11.04 


14:31:57 


0530 


6 Hebe 


319.527 


-23.694 


23.311 


2006.11.07 


06:09:23 


0540 


1 Ceres 


322.469 


-26.277 


160.33 


2006.11.09 


01:08:33 


0560 


93 Minerva 


152.718 


19.193 


4.619 


2006.11.20 


09:13:18 



)i, cr= 26570, 0.0011, 0.0306 



by the two solid lines in Fig.|3}, the mean standard deviation 
(cr) is reduced to 0.025 ± 0.001 ADU. 




Figure 2. The final data histogram for a typical TSD file. # is 
the number of data points, and a is the mean and sigma of the 
fitted Gaussian function respectively. 




200 
TSD file index 



Figure 3. The standard deviation distribution in the AKARI PV 
phase data near the NEP. The median value of cr is 0.025 ADU. 



line in each detector pixel. The smoothing length is equiva- 
lent to an angular size of ~15 arcseconds. 

(4) We estimate the standard deviation and then ex- 
clude any signal with S/N > 3. 

(5) We repeat step (3) until the standard deviation con- 
verges. The final data histogram is fitted with a Gaussian 
curve to derive the mean ^ and the standard deviation o. 

Fig. |2] shows the final data histogram of a typical TSD 
file or scan and its Gaussian fit. Averaged over 412 scans, (tj) 
is 0.026±0.003 ADU. After excluding outliers defined as |cr- 
((t) I > 0.003 ADU (plus signs outside the region delineated 



2.2 90 micron flux interpolation 

The AKARI WIDE-S band channel has a 90 fim effective 
wavelength. To estimate 90 /im fiux, we need to interpolate 
using IRAS fiux measurements at 60 and 100 /^m (Fl^^^ 
and Fwo^^)- According to the IRAS 100 /im fiux quality, we 
divide our objects, all of which have high quality FIq'^^ , 
into two categories: (I) moderate quality or (II) upper limit 
(i.e. undetected). One caveat to bear in mind is that IRAS 
fiux densities with moderate quality are not highly reliable 
as they do not satisfy the reliability requirement of the FSC 
which is ^ 90% at 12, 25 /^m and ^ 80% at 60 yLxn. Moreover, 



moderate quality F-^ 



IRAS 



may be caused by cirrus and not 



the source itself (Moshir et al. 1992). 

Consequently, we find 23 IRAS sources in Cat. (I). The 
interpolated 90 /xm fiux Fgo* is obtained by fitting FIq-^^ 
and FIix)"^^ with au x Bi, (T) curve based on the far-infrared 
spectrum derived from radiative transfer models for cirrus 
and starburst galaxies. In Tabled we present source name, 
non-colour corrected F/o^^^, Fl^i''^ and Fg'J^*, redshift (if 
available) and observed wavebands for each source. The av- 
erage colour at 90-60 iim (Fgo/Feo) is 2.48. Colour correction 
is needed for SEDs other than F^ (x -. For ordinary stars 
and galaxies, this correction is of order a few percent or even 
less. Therefore, due to the relatively large uncertainty in the 
fiux interpolation and the peak timeline signal as a fiux in- 
dicator for point sources, colour correction is neglected in 
this paper. 

To further check the reliability of the IRAS sources in 
Cat. (I), we have used the optical counterpart identification 
method. Fig.|4]shows postage stamps from the Digitized Sky 
Survey (DSS) centred at the IRAS position of each source. 
About 6 IRAS sources (shown in yellow colour) do not have 
obvious counterparts present on the optical photographs. 
From Table |3] we find that these sources (with * in front of 
the source name) have only been observed at 60 and 100 /im 
and their 60 /im fiux densities are ~ 0.2 Jy, i.e. they are very 
faint objects. However, all IRAS sources in Cat. (I) appear to 
be reliable extractions judging from the IRAS 60 and 100 /im 
maps (Fig. [5] and Fig. |6j. The IRAS source F18001-f-6636 
seems to be missing in these maps because of the adjacent 
Cat's Eye Nebula. Fig.[7]shows the co-added AKARI 90 fim 
images centred at the IRAS positions of each source in Cat. 
(I). In Section [3.21 we will derive the conversion factor by 
comparing Fgo' with the AKARI peak timeline signal for 
each IRAS source in Cat. (I). 
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Figure 4. DSS images (image size=0.5 X 0.5 arcmin^; image scale=linear) of the IRAS sources in Cat. (I), centred at the IRAS position 
of each source. The images are ordered in the same way as in Table|3]from left to right, top to bottom. 




Figure 5. IRAS 60 /.tm images (image size=0.5 X 0.5 deg.^; image scale=linear) of the IRAS sources in Cat. (I), centred at the IRAS 
position of each source and ordered in the same way as in Fig. |4] 



For IRAS sources in Cat. (II), i.e. sources undetected 
at 100 /^m, we can predict 90 yLxn flux density by adopting 
(Fioo/Jiso) = 2.24, which is derived from the IRAS Faint 
Source Survey (FSS) redshift survey of 700 square degrees. 
Although IRAS sources in Cat. (II) will not used in calibra- 
tion, it is interesting to compare the predicted flux densities 
with the AKARI measurements. 



3 TIMELINE ANALYSIS 
3.1 Procedure outline 

Timeline Analysis, developed based on the IRAS PSC 
model, is a one-dimensional source detection method. It de- 
tects and outputs SCONs in each scan after deletion of sig- 
nals below a certain threshold. Although the AKARI All- 
Sky Survey has a dedicated point source extraction and pho- 
tometry software package, SUSSEXtractor (Savage & Oliver 
2007), Timeline Analysis is a valuable tool at the initial per- 
formance evaluation stage because its transparency provides 
insight on the timelines and guidelines on parameter tuning 
in the dedicated source extraction pipeline. 

Our goal is to derive the conversion factor which trans- 
forms the AKARI peak timeline signal to the interpolated 
monochromatic 90 pm flux density of a source. Then we can 
estimate the point source detection limit by converting the 



standard deviation in ADU to noise level in Jansky. The 
following is an outline of the Timeline Analysis method. 

Step 1: We prepare a source- list with positions and in- 
terpolated 90 ^m flux densities in the surveyed region, which 
will be used to cross-match with timeline detections later. 

Step 2: We read in all TSD flies in the data set and 
outputs SCONs above a certain threshold, e.g. S/N > 3 or 
5. 

Step 3: SCONs detected in each TSD flle is cross- 
matched with the source-list within a searching radius of 
30 arcseconds, in accordance with the positional accuracy of 
the GADS. The mean RA and DEC value of the matched 
SCONs for a single target source as the AKARI position 
and the peak timeline signal is an indicator of the AKARI 
90 /im flux density. 

Step 4: Finally, the conversion factor is computed as the 
ratio of the interpolated 90 /^m flux F90 * to the peak timeline 
signal for every seconds-conflrmed detection of a particular 
source. Thus, IRAS sources covered by more TSD flies will 
have more derivations of the conversion factor. The average 
conversion factor will be used to determine the noise level 
of the AKARI All-Sky Survey. 
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Figure 6. IRAS 100 /xm images (image size=0.5 X 0.5 deg.'^; image scale=linear) of the IRAS sources in Cat. (I), centred at the IRAS 
position of each source and ordered in the same way as in Fig. |4] 




Figure 7. AKARI 90 /^m images (image size=0.1 X 0.1 deg.-^; image scale=squared) of the IRAS sources in Cat. (I), centred at the IRAS 
position of each source and ordered in the same way as in Fig. |4] 



Table 3. A total of 23 IRAS sources in Cat. (I). 

Name F^^^^^Jy) F^^'iJy) F/«-4S (Jy) 



Redshift 



Observed Wavebands 



F18001+6636 2.04E+00 2.31E+00 2.17E+00 0.026 B, 2MASS (J, ii'i,),IRAS(12,25,60,100 /.tm), ... 

*F18112+6503 2.10E-01 5.94E-01 6.65E-01 - IRAS (60,100 ^lm) 

F18130+6455 2.43E-01 6.34E-01 6.99E-01 - 2MASS (J, i^T^), IRAS (12,60,100 ^tm) 

F18185+6341 2.04E-01 6.53E-01 7.49E-01 - 2MASS (J, ftTj,), IRAS (60,100 Atm) 

F18197+6339 6.01E-01 1.02E+00 1.04E+00 0.027 2MASS (J, ii's),IR.AS(12, 25,60,100 /im), 170 /^m, ... 

F18252+6315 3.89E-01 9.14E-01 9.87E-01 0.084 2MASS (J, ii's), IRAS (25,60,100 Atm) 

F18286+6309 7.72E-01 8.33E-01 7.74E-01 - IRAS (25,60,100 fim) 

*F18293+6304 2.05E-01 5.93E-01 6.67E-01 - IRAS (60,100 ^tm) 

F18223+6255 1.59E-01 5.42E-01 6.30E-01 - 2MASS (J, /sTs), IRAS (25,60,100 ^tm) 

*F18254+6200 1.72E-01 5.95E-01 6.93E-01 - IRAS (60,100 ^J.m) 

F18301+6138 2.40E-01 6.45E-01 7.15E-01 - 2MASS (J, ftTs), IRAS (60,100 Atm) 

F18425+6036 l.OlE+01 1.93E+01 2.01E+01 0.013 2MASS ( J, H, _fs:s),IRAS(12,25,60,100 ^tm), ... 

F18353+5950 5.74E-01 1.17E+00 1.24E+00 0.028 2MASS (J, H, i<'s),IRAS(25,60,100 /xm), ... 

*F18436+5931 2.11E-01 7.80E-01 9.21E-01 - IRAS (60,100 ^J.m) 

F18440+5900 3.37E-01 8.79E-01 9.69E-01 - 2MASS (J, ff, Xs), IRAS (25,60,100 ^tm) 

F18436+5847 5.69E-01 1.04E+00 1.08E+00 - 2MASS (J, Jf, A's), IRAS (25,60,100 ^tm) 

F18506+5801 3.79E-01 9.12E-01 9.89E-01 0.029 B, 2MASS ( J, /«:s),IRAS(25,60,100 /^m), 170 /^m, ... 

*F18520+5715 2.11E-01 7.27E-01 8.46E-01 - IRAS (60,100 ^J.m) 

F19026+5654 2.99E-01 7.23E-01 7.85E-01 - 2MASS (J, ii's), IRAS (25,60,100 a*™) 

*F18597+5631 1.64E-01 5.43E-01 6.27E-01 - IRAS (60,100 Atm) 

F18499+5542 4.28E-01 l.OOE+00 1.08E+00 - IRAS (60,100 ^im), ... 

F18510+5539 1.25E+00 2.41E+00 2.51E+00 0.0242 MASS (J, /fs), IRAS (25,60,100 /^m), 170 /^m, ... 

F19009+5507 5.02E-01 8.28E-01 8.35E-01 - 2MASS (J, ii's), IRAS (25,60,100 Atm) 
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Figure 8. 12 seconds- and hours-confirmed timeline detections of the IRAS FSC sources in Cat. (I), thresholding at 3fr (open circles) or 
5(T (plus signs). The 23 Cat. (I) sources are colour-coded according to their 60 /^m fiux densities (green: FgO > 0.6 Jy; red: 0.3 < Fgo < 0.6 
Jy; blue: Feo < 0.3 Jy). 



3.2 Timeline source detection 

Thresholding at Sct, we have detected 21 IRAS sources in 
Cat. (I) using the Timeline Analysis method and the two 
missing objects are F18293-F6304 and F18254-H6200. To 
increase the reliability of the extracted sources, we have 
discarded 8 sources which are only present in one TSD 
and 1 source whose peak timeline signals in 2 TSDs 
vary by more than a factor of 2. Therefore, the remaining 
12 sources are F18001+6636, F18112-K6503, F18130-H6455, 
F18197+6339, F18252-I-6315, F18286+6309, F18425+6036, 
F18353+5950, F18436+5847, F18506+5801, F19026+5654 
and F19009+5507. In Fig.[Hl we have plotted the 12 seconds- 
and hours-confirmed timeline detections. The IRAS FSC 
sources in Cat. (I) are colour-coded by the 60 micron flux 
measured by IRAS, F^-^^ . At S/N > 3, nearly all sources 
with fIq'^^ > 0.3 Jy (green and red dots with open circles) 
as well as 2 of the faint ones (blue dots with open circles)are 
detected with seconds- and hours-confirmation. At S/N > 5, 
3 out of 5 sources with Fl^'^^ > 0.6 (green dots with plus 
signs in the middle) are seen. Therefore, one can guess the 
equivalent 5a point source detection limit at 60 fim for a 
single detector row is ~0.6 Jy. As the AKARI has covered 
almost the entire sky at least twice and each scan has 3 de- 
tector rows, ideally the 5a point source flux detection limit 
for a single scan will reach ~0.2 Jy which is deeper than 
the detection limit of the IRAS PSC and comparable to the 
FSC. 

Conversion factors derived from the 12 conflrmed 
sources above 3a and asteroids are shown in Fig. [S] For 
sources with FgJ* < 0.7 Jy (left of the dotted line in Fig.E]), 
the conversion factor appears to be systematically lower 
than that from brighter sources. Since IRAS does not pro- 
vide accurate flux determinations for these faint objects, we 



An investigation of the timelines shows that these sources are 
usually covered by only one TSD or at the end points of the TSD . 



will exclude them in estimating the mean conversion factor. 
Source F18425+6036 is a nearby extended object and hence 
is excluded as well. So, we are left with a total of 9 Cat. 
(I) IRAS sources used for absolute calibration. The mean 
conversion factor is computed as a weighted average, 



(1) 



ai is the uncertainty in the conversion factor, taking into 
account the scatter in Jgo* and the AKARI peak timeline 
signal. The uncertainty in the estimate of the conversion 
factor is calculated as 



^{CF) = 1/ 



(2) 



Therefore, combined with the robust noise estimation in Sec- 
tion[2Tl we deduce that in the AKARI WIDE-S 90 fj.m band, 
the point source flux detection limit for a single detector row 
at S/N > 5 is 1.1 ± 0.1 Jy (or ~0.7 Jy at S/N > 3). 

Fig. [To] shows a comparison between the AKARI 90 ^m 
flux, i.e. peak timeline signal x conversion factor, and the 
interpolated 90 fim flux for Cat. (II) sources detected at 
a 3a threshold. For sources with Kjy* less than ~0.7 Jy, 
the consistency between Fgp ' and Fg^^"^^' begins to break 
down. This trend seems to coincide with the boosting of the 
AKARI fluxes at the faint end present in Fig. [9] It could 
be caused by the well-known flux overestimation effect from 
upward noise fluctuations (Moshir et al. 1992), non- linear 
absolute calibration, incorrect IRAS fluxes and flux interpo- 
lation (Jeong et al. 2007), or by mistaking correlated noise 
for faint point sources as we are approaching the detection 
limit of ~0.7 Jy at S/N > 3. 
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Table 4. 6 cirrus candidates detected at ia or 5o". 
NAME RA (J2000) DEC (J2000) 



CI 


277.796 


61.716 




C2 


286.906 


56.427 


ctor 


C3 


284.842 


57.186 


C4 


276.261 


63.704 


D 


C5 


286.365 


55.464 


C 

o 


C6 


286.332 


55.483 





3.3 Non-IRAS sources and cirrus emission 

The expected number density of sources at ~0.5 Jy is 1 
source per square degree from 90 /xm integral source count 
models based on the European Large Area ISO Survey (Ef- 
stathiou et al. 2000; Heraudeau et al. 2004). In our data set, 
each scan is about 3.6 degree long and 8.2 arcmin wide. In 
other words, each scan covers an area of ~0.5 square degrees. 
At a 5(T threshold, two overlapping scans usually give rise 
to 1 - 2 seconds- and hours-confirmed sources. Therefore, 
with the current high source number density and sensitivity 
estimate, a large fraction of the detections are expected to 
be spurious point sources. In principle, false detections can 
be caused by cosmic ray glitches, transient behaviour of the 
detectors, correlated noise, or filamentary structures such as 
interstellar dust clouds or cirrus discovered by IRAS (Low 
et al. 1984) which dominate the background radiation in the 
far-infrared. Jeong et al. (2005) estimated the sky confusion 
noise (Helou & Beichman 1990; Gautier et al. 1992) due to 
cirrus emission using simulated high-resolution dust maps 
for space missions such as ISO, Spitzer, AKARI, Herschel 
etc. Jeong et al. (2006) included source confusion in different 
source distribution models in determining the far-infrared 
point source detection limit. 

To investigate the effect of infrared cirrus emission on 
source extraction, firstly we need to select seconds- and 
hours-confirmed detections with a high signal-to-noise ratio 
to decrease the contamination from transient signals or mov- 
ing objects. Having selected 5 pairs of overlapping scans, we 
output SCONs above a certain threshold (varies between 4(t 
and 5cr) and then discard those without hours-confirmation. 
Thus, we have obtained a sample of 8 seconds- and hours- 
confirmed sources which are not IRAS sources. In Fig. 1111 
we have overlaid 6 of them (white open circles) on the IRAS 
100 /^m emission maps (Schlegel, Finkbeiner & Davis 1998) 
and found that they generally reside in high or intermediate 
intensity environments. Fig. [12] shows the co-added AKARI 
90 images with greater resolution and confirms the ex- 
tended structures of the selected non-IRAS sources shown 
as green open circles (could be compared with the 90 /xm im- 
ages of the point sources in Fig. [7]l . From now on, we refer to 
these objects as cirrus candidates. In Table |4l we list name 
and position for each candidate. In Section [4.51 we will use 
wavelet multiscale decomposition technique to verify or fal- 
sify the assumption that these seconds- and hours-confirmed 
non-IRAS sources are due to cirrus emission. 

For a fair comparison with the cirrus candidates selected 
at S/N > 4 or 5, we have run our timeline analysis method 
with a threshold of Ag. The extracted seconds- and hours- 
confirmed IRAS sources at S/N > 4 are listed in Table [5] 
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Figure 9. Conversion factor (Fgg'^ / peak timeline signal) vs. 
the interpolated 90 fim flux -FgQ ' (asterisks: IRAS FSC sources in 
Cat. (I) detected at S/N > 3; triangles: asteroids). The dotted 
line is where the interpolated 90 fim flux is equal to 0.7 Jy. For 
clarity, error bars are not plotted. 
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Figure 10. The interpolated 90 /xm flux vs. the AKRAI 90 /^m 
flux for IRAS FSC sources in Cat. (II) detected at S/N > 3 (plus 
signs) and asteroids (triangles). 



4 WAVELET TRANSFORM AND SOURCE 
DETECTION 

4.1 Key motivations 

Infrared astronomical data are often complex combinations 
of random noise (photon shot noise, instrumental noise, 
etc.), variable sky background radiation, point and quasi- 
point sources (e.g. asteroids, stars, double stars, galaxies, 
...), different types of glitches caused by cosmic rays, ex- 



Table 5. 3 IRAS sources in Cat. (I) detected at 4cr. 

NAME RA (J2000) DEC (J2000) 

F18001+6636 270.042 66.6118 
F18197+6339 275.029 63.679 
F18436+5847 281.115 58.837 
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Figure 11. IRAS 100 /im maps (image size=2 X 2 deg.'^; image scale=histogram). From left to right, the white open circles (diameter 
= 4.8 arcminutes, close to the IRAS resolution) are the cirrus candidates CI, C2, C3, C4, C5 and C6 respectively (C5 and C6 are shown 
together in the rightmost map.). 




Figure 12. AKARI 90 /xm maps (image size=0.2 X 0.2 deg.'^; image scale=histogram). From left to right, green open circles (with 
diameter = 1.2 arcminutes) are the cirrus candidates CI, C2, C3, C4, C5 and C6 respectively (C5 and C6 are shown together in the 
rightmost map.). 



tended structures such as galaxies, clusters, cirrus clouds, 
etc. Different signals prefer to exhibit themselves on differ- 
ent scales because they originate from a hierarchy of phys- 
ical structures. For instance, noise tends to dominate on 
small scales, infrared cirrus emission often shows up on large 
scales, while stars and galaxies would appear on somewhat 
intermediate scales. So, it is advantageous to use multires- 
olution techniques to detect signals from objects with a 
particular size. In addition, because our data always con- 
tain some discontinuities and sharp features, the traditional 
Fourier transform becomes an inadequate and inefficient 
analysis tool for modelling the real data. 

For the AKARI All-Sky Survey, a multiscale analysis 
can effectively eliminate spurious sources if the majorities 
have a different spatial scale from the point source scale. 
Essentially, we use wavelet multiresolution decomposition as 
a means of understanding data characteristics and increas- 
ing the reliability and completeness of the future catalogues 
from the AKARI All-Sky Survey. 



4.2 Introduction to wavelets and wavelet 
transforms 

Wavelet transform is in many ways similar to Fourier trans- 
form. Both are linear operators which decompose signals 
onto a set of basis functions. The fundamental difference 
is that wavelet transform adopts localised basis functions 
called wavelets instead of infinitely ranging sinusoidal func- 
tions which have no positional resolution at all. The mother 
wavelet function gives rise to a family of wavelet functions 
by operations of translation and scaling without changing its 
shape. Therefore, wavelet transform offers a time-scale two- 
dimensional representation of the data, i.e. an accurate lo- 
cal description of the frequency components present in each 



signal segment. It is particularly useful for analysing non- 
stationary signals whose characteristics change over time. 
Alternatively, it can be viewed as a multiresolution decom- 
position which is perhaps more widely used (refer to Burrus, 
Gopinath & Guo 1998 for a mathematical background of 
wavelet analysis and Starck & Murtagh 2002, van der Berg 
2004 and references therein for applications of wavelets in 
astronomy and physics). 

There are two types of wavelet transform, the contin- 
uous wavelet transform (CWT) and the discrete wavelet 
transform (DWT). Theoretically the CWT is infinitely 
redundant and inefficient, however it provides a smooth 
wavelet transform of the signal and is suited for feature 
recognition. The DWT is fast to compute and the signal 
reconstruction is straightforward. In this paper, we use the 
CWT to analyse the TSDs. 

The CWT of a square-integrable function f{x) is de- 
fined as 

W{a.,b) = ^ /(a,)^*(^)da:, (3) 

where ^*{x) denotes the complex conjugate of tp{x). By the 
convolution theorem, the wavelet transform is the inverse 
Fourier transform of the product of f{w) and 7/)*(aw), 

/ + 00 
^'■{aw)f{w)e''""dw. (4) 
- oo 

The wavelet power spectrum (WPS) is defined as |M^(a, &)|^. 
Regions of large or small power in the WPS plot indicate 
important or negligible features within the signal. 

There are several factors in determining which wavelet 
function to use (Torrence & Compo 1998). Essentially, 
wavelet transform calculates the correlation of the wavelet 
function and the local signal. Therefore, to detect point 
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sources in the timelines, we need wavelet functions which 
can produce large wavelet coefficients at the position of the 
source and at the same time reduce the effect of noise. We 
use the Gaussian second derivative (also known as the Marr 
wavelet or Mexican Hat wavelet) as our analysing wavelet 
function. It has been suggested by many that the optimal 
wavelet to detect point sources is perhaps the isotropic Mex- 
ican Hat wavelet (Cayon et al. 2000; Gonzalez-Nuevo et al. 
2006). 



coefficient to zero if Wa.b < ka{a) (insignificant wavelet co- 
efficient) . k is chosen to be 2 in this paper. This is known as 
hard thresholding. 

5. Finally we use seconds- and hours-confirmation 
to discriminate real sources from spurious ones. This is 
achieved by comparing the wavelet transform after thresh- 
olding in different detector rows (seconds-confirmation) and 
different TSDs (hours-confirmation). 



4.3 Wavelet multiscale representation: definition 
and outline 

On what scales do we expect point sources to appear? In the 
WIDE-S 90 fim band, the FWHM of the PSF is ~39 arcsec- 
onds, the sampling rate is 25.28 Hz, and the scan speed is 
~3.6 arcmin per second. Therefore, a single point source can 
be observed by ~9 samplings in one detector pixel provided 
no data points are fiagged out, 

..Im'"'^'''/ >^ 25.28 H. . 9. (5) 
3.6 X 60 arcsecf sec 

Therefore, the closest wavelet width to 9 is 2'^ = 8 under a 
dyadic sampling strategy. In other words, the optimal res- 
olution leveQ for detecting point-like objects is a = 2^. It 
follows that quasi-point sources or slightly extended sources 
mostly inhabit on scale 2*. To sum up, our prior knowledge 
of point source detection in the multiscale representation 
model includes: 

PI. At resolutions 2^ and 2^, signals are assumed to be 
dominated by random noise; 

P2. Point sources can be most easily detected at res- 
olution 2'^. Quasi-point sources / slightly extended sources 
prefer to show up at resolution 2*. We will refer to these 
resolution scales (2^ and 2^) as the Point Source Scales; 

P3. At resolutions coarser than 2*, cirrus emission and 
the sky background begin to dominate. 

We use the IRAS sources in Table [5] to demonstrate 
source detection with wavelet multiscale decomposition 
technique and also to compare with the cirrus candidates. 
A summary of the procedures is the following. 

1. In each detector row[f|, we search for the detector pixel 
which passes a target source. The 3 passing detector pixels 
should be approximately 20 or 21 detectors away from each 
other due to the configuration of the detector arrays of FIS 
and the rotation angle (26.5 deg) with respect to the scan- 
ning direction. 

2. In each detector pixel, the number of data points A'' is 
required to be a power of 2, i.e. N ^2-' . Therefore, we add 
an adequate amount of zeros at the end of each timeline. 

3. We compute the CWT of the timelines on the Point 
Source Scales (2^ to 2'*), producing a set of wavelet coeffi- 
cients 'Wa,b on each scale. 

4. We calculate the standard deviation o-(a) of the 
wavelet coefficients at each scale and then set any wavelet 



* This is the Object Scale defined as the scale where the maxi- 
mum wavelet coefficient of an object lives in the Multiscale Vision 
Model (Bijaoui & Rue 1995; Stark & Murtagh 2002) 
^ The FIS has 3 detector rows in the WIDE-S band, each of which 
has 20 detector pixels (Table 1). 



4.4 Multiscale analysis of the IRAS sources 

Examples of the original timelines in one scan for the 3 IRAS 
sources in Table 5 are shown in Fig. [T31 For each source, 
there are 3 timelines in the 3 passing detector pixels respec- 
tively. The position of the target source is indicated by the 
intersection between the signal and the vertical dotted line. 
The brightest source F18001-I-6636 withKjj* = 2.3 Jy can 
just be detected by eye, while the other two sources with 
^90* = 1-0 Jy are obscured by the noise. 

In the same figure, we have shown the continuous 
wavelet transform of the timelines after thresholding for each 
IRAS source in 2 overlapping TSDs (shown as TSD 1 and 
TSD 2) on scales 2^ and 2'*. In fact, the reconstructed time- 
line on each scale is proportional to the amplitude of the 
wavelet coefficient (Torrence & Compo 1998) . The timelines 
in the 3 passing detector pixels are represented by 3 differ- 
ent colors, black, red and blue. If more than two positive 
signals of different colours overlap at the same position in 
one TSD, then it is regarded as a seconds-confirmed detec- 
tion. To discriminate real sources from spurious ones, we also 
require hours-confirmation by checking the other TSD at the 
same position. Only when a seconds-confirmed detection is 
present in both TSDs, i.e. a seconds- and hours-confirmed 
source, do we acknowledge it as a real source. 

So, the first conclusion from Fig. [13] is that each IRAS 
source can be clearly identified as the only seconds- and 
hours-confirmed detection in the entire TSD at resolution 
2'*. Secondly, although F1800H-6636 is also present with 
seconds- and hours-confirmation at resolution 2^, the two 
fainter sources F18197-I-6339 and F18436+5847 are not de- 
tected according to our selection criteria at j = 3. In addi- 
tion, due to the noisy nature of the data, F18197-I-6339 and 
F18436+5847 will be missing at j = 4 if a higher threshold 
is used. 

Fig. [2] shows the WPS as a function of position and 
scale for IRAS source F18197-I-6339. Again, there are 6 pan- 
els (3 passing detector pixels per TSD x 2 TSDs) and the 
vertical dotted line in the middle of each panel marks the 
position of the target source. In order to compare the WPS 
in different detector pixels and for different objects, we use 
common contour levels in Fig. 1141 Fig. 1171 Fig. [18] and Fig. 
1191 The wavelet power increases towards red in each contour 
plot. A high power at the position of the IRAS source and 



on the correct scales (around 2 



can be unmistakably 



identified in almost every contour panel of F18197-I-6338. 
The missing power in the rightmost panel in TSD 2 is due 
to the fiagged data around the source (see the top left panel 
in Fig. I13|l . The noisy nature of the data can been seen in 
the contour plots of the wavelet power. For example, there 
are contours of large power in the timelines on a range of 
scales including the Point Source Scales. 
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Figure 13. Left: examples of the original timelines in the 3 passing detector pixels in one TSD for each IRAS source. Right: the 
continuous wavelet transform of the entire timeline in each passing detector pixel after thresholding on the Point Source Scales in 2 
overlapping TSDs for each source. The dotted line shows the position of the target source. 



4.5 Multiscale analysis of the cirrus candidates 

The same wavelet transform procedure has been applied to 
the 6 cirrus candidates selected with S/N> 4 and the re- 
sult is shown in Fig. 1151 Apart from candidate C2, all the 
other cirrus candidates are clearly not seconds- and hours- 
confirmed detections on the Point Source Scales. Therefore, 
according to our selection criteria these sources are spurious 
and not new discoveries by AKARI. 

In Fig. 1161 the original timelines of cirrus candidate C2 
in 2 TSDs (TSD 1 and TSD 2) are shown. The bottom and 
top timeline in TSD 1 and the middle timeline in TSD 2 
clearly show the presence of a spike at the position of the 
cirrus candidate. The contour plots of the WPS of C2 are 
shown in the lower half of Fig. [T7] In both TSDs, the effect 
of a spike is manifested by contours of power resembling 
that of point sources at the vertical dotted line (the left and 



right panel in C2 TSD 1 and the middle panel in C2 TSD 
2). For other cirrus candidates, contour plot similar to that 
of a point source can be seen in some detector pixels. An 
investigation of the original timelines shows that it is also 
due to spikes found at the position of the cirrus candidates. 

If these non-IRAS detections detected in Timeline Anal- 
ysis with seconds- and hours-confirmation are induced by 
extended structures such as cirrus clouds, they should ap- 
pear as confirmed detections on scales larger than the Point 
Source Scales. However, for every cirrus candidate studied 
in this paper, this assumption can be rejected based on the 
WPS plots from Fig. [T7] to Fig. [T^ In some contour plots 
a thin curved line can be seen at either side of the panel 
(CI, C2 and C3). Regions above the thin line are the cone 
of influence (Torrence & Compo) where the edge effect of 
the TSD makes the power spectrum unreliable. In the case 
of cirrus candidate CI, regions of high power on large scales 
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Figure 14. The local wavelet power spectrum. The x-axis is centred at the position of the target source (the vertical dotted line). 
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Figure 15. The continuous wavelet transform after thresholding for cirrus candidates CI to C6. Each panels shows the wavelet coefficients 
in 3 passing detector pixels in 2 overlapping TSDs. The position of the cirrus candidate is indicated by the dotted line. 
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Figure 16. The original timelines of cirrus candidate C2 in the 3 passing detector pixels in 2 overlapping TSDs. 



are present in the left panel of TSD 1 as well as TSD 2 at 
the vertical dotted line, while they are absent in other de- 
tector pixels. Similarly, for cirrus candidates C4, C5 and C6, 
at the position of the candidate, only one detector pixel dis- 
play high power on large scales in each case (the right panel 
in C4 TSD 2; the left panel in C5 TSD 2; the middle panel 
in C6 TSD 2). For C2 and C3, apart from spikes seen in 
some detector pixels, the dotted linen in does not intersect 
with regions of large power. The fact that the contour plots 
of wavelet power along the dotted line look wildly different 
across different detector pixels of the same object leads us 
to conclude that these spurious sources are not caused by 
cirrus emission but random noise in the timelines. 



5 DISCUSSIONS AND CONCLUSIONS 

The AKARI All-Sky Survey has covered over 90 per cent 
of the entire sky twice. So, theoretically the 5a point source 
flux detection limit at the WIDE-S 90 //m band will reach 1.1 
Jy / -^2 {overlapping scans) x 3 {detector rows) w 0.4 
Jy assuming equal data quality among different detector 
rows and different scans, which will give rise to more than 
40,000 extragalactic sources over the whole sky. In regions 
covered by multiple scans (e.g. the NEP), the point source 
detection limit will be improved further. Therefore, future 
catalogues from the AKARI All-Sky Survey will be valuable 
samples for large-scale structure studies. Further improve- 
ment on the absolute calibration can be made with better 
quality data from the data reduction pipeline which is under 
intensive development, larger data sets and well- modelled 
objects (comets, stars, galaxies etc.). 

The localising property of wavelet transform prevents 
transient signals being submerged in the average character- 
istics of the overall timeline. Wavelet transform, as a band- 
pass filter, decomposes data onto a succession of resolution 
levels and hence separates signal components due to phys- 
ical structures of different sizes. In our wavelet multiscale 
analysis, IRAS sources are the only detections with seconds- 
and hours-confirmation at the expected resolutions for point 
sources or small extended sources. However, as our data set 
only provides us with 3 sources above a 4a threshold, a larger 
sample will be desirable to further test our method. 

By thresholding at a high signal-to- noise ratio (S/N > 4 



or 5), we find that most of the seconds- and hours-confirmed 
sources which are not in the IRAS point source catalogues 
trace cirrus clouds, i.e. they reside in environments with 
strong or intermediate 100 /um emission. In the co-added 
AKARI 90 ^m maps, it is clearly shown that these non-IRAS 
sources referred to as cirrus candidates are associated with 
bright extended structures, either in the in-scan direction 
(the time direction) or both the in-scan and cross-scan direc- 
tion. Unlike the IRAS sources, most of the cirrus candidates 
fail to to be detected with seconds- and hours-confirmation 
on the Point Source Scales under continuous wavelet trans- 
form. It confirms that they have a different scale preference 
compared with the IRAS sources. However, these non-IRAS 
point sources are not likely due to cirrus judging from the 
contour plots of the wavelet power. At the position of each 
cirrus candidate, usually one or two detector pixels exhibit 
high power on large scales, while other detectors either show 
presence of spikes on the Point Source Scales or no regions 
of high power at all. In other words, the WPS differ signifi- 
cantly across different detector pixels in the same scan and 
also different scans of the same candidate. With improve- 
ment of data processing currently in progress, we can filter 
out excessive noise in the raw timelines to increase the relia- 
bility and completeness of the AKARI All-Sky Survey data 
products. 
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